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FOREWORD 


This report covers the work completed on the research project "Effects 
of Angle of Attack on the Coupled Radiative and Convective Heat Transfer 
About Blunt Planetary Entry Bodies" during the period May 1 to October 31, 
1978. The work was supported by the NASA/Langley Research Center (Aero- 
thermodynamics Branch of the Space Systems Division) through research grant 
NSG 1464. The grant was monitored by Dr. Randolph A. Graves, Jr. of the 
Space Systems Division. 
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LAMINAR AND TURBULENT FLOWS OVER SPHERICALLY BLUNTED 
CONE AND HYPERBOLOID WITH MASSIVE SURFACE BLOWING 

By 

Ajay Kumar 1 and S. N. Tlwarl 2 
SUMMARY 

Numerical solutions are presented for the flow over a spherically 
blunted cone and hyperboloid with massive surface blowing. Time-dependent 
viscous shock-layer equations are used to describe the flow field. The 
boundary conditions on the body surface include a prescribed blowing-rate 
distribution. The governing equations are solved by a time-asymptotic 
finite-difference method. Results presented here are only for a perfect 
gas-type flow at zero angle of attack. Both laminar and turbulent flow 
solutions are obtained. It is found that the effect of the surface 
blowing on the laminar flow field is to smooth out the curvature discon- 
tinuity at the sphere-cone juncture point, which results in a positive 
pressure gradient over the body. The shock slope increases on the 
downstream portion of the body as the surface blowing rate is increased. 

The turbulent flow with surface blowing is found to redevelop a boundary- 
layer-like region near the surface. The effects of this boundary-layer 
region on the flow field and heating rates are discussed. 

1 . INTRODUCTION 

The exploration of the outer planets requires the development of a 
reliable aerothermal environment to be encountered by the entry probe. 

The aerothermal environment is characterized by high energy flow, large 
heat-transfer rates to the probe's surface, and high rates of mass injection 
from the probe's ablative heat shield into the flow. This requires the 
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development of computer codes which can provide coupled solutions for chemi- 
cally reacting and radiating flow with massive ablating rates. Most of the 
coupled solutions presented so far are based on the flow over analytical 
bodies such as hyperboloids (ref. 1 to 3). These computer codes have 
difficulty In treating the realistic entry probe shape, viz, the spherically 
blunted cone with surface blowing. Reference 4 presents some results with 
surface blowing on spherically blunted cones for air, but the blowing rates 
considered are very low, and the blowing-rate distribution Is such that It 
starts at the stagnation point and goes to zero at the sphere-cone juncture 
point. There Is no blowing on the conical portion of the body. The purpose 
of this paper Is to present the solutions for the flow over a hyperboloid 
and spherically blunted cone with massive surface blowing rates. Results 
presented here are for a perfect gas-type flow at zero angle of attack and 
do not Include radiative heating or chemical factions in the flow field. 

The same gas Is injected from the surface as is In the main flow field. 

The analysis of reference 5 is modified and usti here. Time-dependent 
viscous shock-layer equations in the body-oriented coordinate system are 
used to describe the flow field. A time-dependent finite-difference method 
due to MacCormack (ref. 6) is used to solve the equations. The boundary 
conditions on the body surface include the surface blowing. A special 
central differencing of the type suggested in reference 7 is used in the 
tangential direction at the sphere-cone junction point to take into account 
the discontinuity in the surface curvature. Shock is treated as a sharp 
discontinuity across which shock relations are used to compute the flow 
quantities behind the shock. The effects of turbulence on the flow-field 
quantities are also studied in the presence of surface blowing. A two- 
layer eddy viscosity model, consisting of an inner law based upon Prandtl's 
mixing length concept and the Clauser-Klebanoff expression for the outer 
law, is used in the analysis to approximate the eddy viscosity. 

It is found that the effect of the surface blowing on the laminar 
flow field over a sphere-cone is to smooth out the curvature discontinuity 
at the juncture point which results in a positive pressure gradient over 
the body. Shock standoff distances also exhibit similar behavior. The 
shock gets steeper and steeper on the downstream portion of thd body as the 
blowing rate is increased. The laminar shock layer over the sphere-cone 
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In the presence of surface blowing Is characterized by an Inner Invlscld 
layer* a thin shear layer, an entropy layer, and an outer Invlscld layer. 

The turbulent flow with surface blowing Is found to redevelop a 
boundary-layer-1 Ike the region near the surface. The shock standoff 
distances are significantly reduced over the downstream part of the body, 
and, unlike laminar flow, the shock does not get steeper with increasing 
blowing. The Increased gradients near the surface In the turbulent flow 
make the convective heating and skin friction significant even with massive 
blowing, whereas these quantities are negligibly small for the corresponding 
laminar flow. Another Important effect of the formation of the boundary 
layer like the region near the surface is on the radiative heat flux toward 
the body. Reference 2 shows that the radiation blocked in the Inner layer 
Is greatly reduced for the turbulent flow, thus significantly Increasing the 
radiative heat transfer to the body surface. 

2. NOMENCLATURE 
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nondimens lonal total enthalpy, H'/IT 2 

nondimensional specific enthalpy, h'/lT 2 

free-stream Mach number 

surface blowing rate, (p * v') w /(p' V) 

surface blowing rate at the axis of symmetry 

coordinate normal to the body, n'/R^ 

nondimensional pressure, p’/p^, IT 2 

free-stream pressure (N/m 2 ) 

gas constant ( J/kg K) 

free-stream Reynolds number, p^ R^/y^ 

nose radius (m) 

body radius normal to the body axis, r ' /R^ 
coordinate measured along the body, s'/R^ 
nondimensional temperature, T'/T ^ 
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free-stream temperature (K) 
nondimensional time, t* V'/Ri 
transformed time 

nondlmens lonal tangential velocity, u'/lT 

nondimensional normal velocity, v7\T 

free-stream velocity (m/sec) 

transformed coordinate given by equation (3) 

tranformed coordinate given by equation '6) 

transformed coordinate given by equation (3) 

transformed coordinate gwen by equation (6) 

r + n cos 9 

stretching factor 

ratio of the specific heats 

shock standoff distance, S'/R^j 

body angle measured from the body axis 

Prandtl number 

turbulent Prandtl number 

1 + n>c 

local curvature, < ' 
nondimensional density, p'/pj, 
free-stream density (kg/m 3 ) 
nondimensional viscosity, u'/u' 

<x 

nondimensional eddy viscosity, i4/u' 

i ® 

free-stream viscosity (N sec/m 2 ) 
normalized eddy viscosity, uj/u 
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Superscript 

(') dimensional quantities 

Subscript 

w conditions at the body surface 

- conditions In the free stream 

3, ANALYSIS 

The basic equations are obtained from the unsteady Navier-Stokes 
equations by keeping terms up to second order in the inverse square root of 
Reynolds number In both the viscous and invlscld regions. The body-oriented 
coordinate system shown in figure 1 is used. These time-dependent viscous 
shock-layer equations, developed in reference 5, are modified to include the 
turbulence and are used here to describe the flow field. The present analysis 
is for a perfect gas flow at zero angle of attack and does not include radia- 
tive heating or chemical reactions in the flow field. 

3.1. Governing Equations 

The governing equations are expressed in the nondimensional conservative 
form as: 
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These equations are not valid at the axis of symmetry. A limiting form 
of the governing equations is obtained by differentiating equation (1) with 
respect to s and taking the limit as s->0. The following equations are 
obtained at the axis of symmetry: 
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In addition to the above equations, an equation of stats and a viscosity 
law are used to complete the set of governing equations. 

Two independent variable transformations are applied to the governing 
equations. The first transformation maps the computational domain into a 
rectangular region in which both the shock and the body are made boundary 
mesh lines. This transformation Is given by 

y * s, z * 1 - n/5 and t * t (3) 

whore 6 is the shock standoff distance. The transformed equations (1) and 
(2) are 
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* y 
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where 




_ 3 _ 

3Z 


[(1 - Z) <« t u 0 + s y M 0 ) - N 0 ] + «q 0 - 0 


. (5) 
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The computational region is further mapped to another plane to allow 
higher resolution near the body surface without too much increase in the 
number of mesh points in the normal direction. The higher resolution near 
the body surface is desirable for an accurate calculation of skin-friction 
and heating rates. The second transformation is given by (ref. 8) 


y = y, z = 
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With this transformation, the 


and t * t 


final form of equation (4) is 
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= 0 


where 

U' = 611, M' = 6M 

N 1 = (—) [(1 - z) (fir U + fi- M) - N] 

\3z / t y 

Q* = 6Q - ( g.r - J- p) CO - z) (fi T U + 6^ M) - N] 
Equation (5) can also be written in a similar form; here 
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and 


2 * 6 ' [(£-H)* - + ’] 

The partial derivative with respect to n in the expressions for N 
and Q are to be replaced by 

2-a.li_._l (*L\ 3 

an <$ az *6 \az/ az 

3.2. Boundary Conditions 

The boundary conditions at the shock are calculated by using the shock 
relations. At the wall, no slip conditions are used. The wall temperature 
is taken as a specified value, and a prescribed blowing rate distribution is 
used on the surface. 

3.3. Eddy Viscosity Approximation 

A two-layer eddy viscosity model consisting of an inner law and an 
outer law is used in the present analysis. This model, introduced by 
Cebeci (ref. 9), assumes that the inner law is applicable for the flow from 
the wall out to the location where the eddy viscosity given by the inner law 
is equal to that of the outer law. The outer law is then assumed to be 
applicable for the rest of the shock layer. The expressions for the inner 
and outer laws are given below. 

3.3.1. Inner Eddy Viscosity Approximation. -Prandtl's mixing-length 
concept is stated in nondimensional variables as 


p«. 2 Re 3u 

(7) 

u 3n 



The mixing length i is evaluated from the expression 
i = k ! n[l - exp (-n + /A + )] 
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where 



Here, ki is the Von Karman constant, which is assumed to have a value of 
0.4, and A + is a damping factor which is expressed as 

A + * 26 exp (-5.9 v + ) 


where 
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v w /u T 


and 



2.3.2. Outer Eddy Viscosity Approximation. - For the outer region of the 
viscous layer, the eddy viscosity is approximated by the Clauser-Klebanoff 
expression (based on refs. 10 and 11): 
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and 


y. * [1 + 5.5 (n/6.) 6 ] 1 
• » n ° 

Here, 6^ is the boundary-layer thickness and u g is the local value at the 
edge of the boundary layer. The expression for y^ n is Cebeci's approxima- 
tion of the error- function definition presented by Klebanoff (ref. 11). The 
boundary-layer thickness 5^ is assumed to be the value of n at the point 
where 


~ = 0.995 


In the present analysis, the turbulent Prandtl number is assumed to be 

0.9. 


4. METHOD OF SOLUTION 

MacCormack's two-step finite-difference method (ref. 6) is used to solve 
the governing equations. This explicit method has second-order accuracy in 
both space and time and is highly efficient on the CDC-STAR-100 computer. 

Since only the steady-state solution is of interest in the present investi- 
gation, the solution is marched in time on each mesh point according to its 
largest possible Courant-Friedrich-Lewy (CFL) time-step size. It is shown 
in reference 12 that the use of local time-step results in a speedup of two or 
or more over the regular MacCormack scheme in which the solution is marched 
with global minimum CFL time-step on all the mesh points. Due to these 
reasons, the present analysis still uses the original MacCormack scheme 
although several more recent versions are available. 

In addition to the instability which will occur if the maximum allowable 
time-step is exceeded, it is found that the present solutions exhibited 
large oscillations across the shock layer. A fourth-order damping technique 
is used to damp these oscillations. The details of the method are given 
in reference 5. 
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A special central differencing of the type suggested In reference 7 Is 
used In the tangential direction at the sphere-cone juncture point to take 
into account the discontinuity In the surface curvature. 

5. DISCUSSION OF RESULTS 


Flow-field results are presented for the forebody of hyperboloids and 
spherically blunted cones with massive surface blowing. Results are presented 
for both laminar and turbulent flow. Two different sets of conditions are 
used. These conditions are given below. 


I. Gas : Air 


M = 10.33 
00 

T' = 46.26 K 

oo 

= 100.77 N/m 2 


T,' = 330.6 K 
w 

y = 1.4 
a * 0.7 


R' = 0.03175 m Re = 115800 

N 

u’ = 1.458 x 10' 6 {(T'^-VOlO + T')} N-sec/m 2 

R * 287 J/k K 
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The surface blowing rate distribution used with these conditions is given by 


m/m 0 * (P v )y/(P V ) W , 


s 


= 0 


2 


0.3 s + 
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These conditions and surface blowing rates are identical to those used in 
reference 3 and are used to obtain the laminar flow-field results over a 45 
degree half-angle hyperboloid only. 

II. Gas : Hydrogen-Helium Mixture (under perfect gas assumption) 

M =43.84 T' = 4000 K 

oo W 

T' = 145 K y = 1.224 

00 

p' = 1.27E - 4 kg/m 3 o = 0.72 

oo 

R,', = 0.222 m Re = 156700 

N 
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N-sec/m 2 


-7 0.65 

w' - 2.48 x 10 7 (T') 

R - 3593.6 J/ k g K 

These conditions correspond to typical Jovian entry conditions. The surface 
blowing rate distribution used with these conditions is shown in figure 2. 

This set of conditions is used with both the spherically blunted cone and 
hyperboloid to obtain the laminar and turbulent flow-field results. 

Results of the first set of conditions are discussed first. Figure 3 
shows the effect of surface blowing on the shock standoff distance. It is 
seen from this figure that the shock standoff distance increases with 
increasing blowing rates and the shock becomes progressively steeper for 
the downstream points. Figure 4 shows the corresponding surface pressures 
which also increase with increasing blowing rates for the downstream points. 
Thus, the net effect of the mass injection is to produce effectively a 
blunter body. 

The results of the present analysis are compared with those of reference 
3. It is seen from figures 3 and 4 that the results of the two analyses agree 
very well for zero blowing but the difference increases with increasing 
blowing rates. For m Q = 0.8, the difference is the order of 30 percent in 
both the shock standoff distances and surface pressures. No reason could 
be attributed to these large differences. Both the analyses use the same 
governing equations and input conditions. The only difference is in the 
method of solution. 

It is mentioned in reference 3 that the viscous shock-layer equations 
are not extendable to massive blowing rates and that it is necessary to use 
full Navier-Stokes equations with massive blowing. It is found in the 
present analysis that the viscous shock-layer equations can very well handle 
the massive surface blowing rates. In fact, the results obtained in the 
present analysis with viscous shock-layer equations and full Navier-Stokes 
equations are almost identical for the entire range of m Q from 0 to 0.8. 

Figures 5 to 7 show the profiles for tangential velocity, enthalpy, 
and normal velocity at s = 0.9. Figures 5 and 6 show the typical behavior 
of the shock layer in the presence of surface blowing. It is characterized 
by an inner inviscid layer, a thin shear layer in which most of the changes 
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occur, and an outer invlscld layer. The extent of the inner invlscld layer 
increases with the increasing blowing rates. In the case of an ablating 
surface, this inner inviscid layer will mainly be containing the ablation 
gas. 

The normal velocity profiles In figure 7 show that the normal velocity 
decreases slightly In the inner Invlscld layer, Increases in the shear layer, 
and then decreases monotonically to the shock value in the outer inviscid 
layer. This behavior of the normal velocity has been reported in reference 3. 
The Increase In the normal velocity In the shear layer would enhance the 
mixing of the Inner and outer layer gases, thus Increasing the thickness of 
the ablation gas layer. This Increase in the ablation-layer thickness is 
particularly important since most of the attenuation of the radiation takes 
place in this region. 

Part of this section will now discuss the results from the second set 
of conditions. Results of the laminar flow are discussed first. Figure 8 
shows the surface pressure distribution over a 45-degree half-angle sphere- 
cone with increasing surface blowing. It is seen from this figure that for 
zero blowing there is a negative pressure gradient on the conical portion 
of the body due to the discontinuity in the curvature at the juncture point. 

As the blowing rate is increased to m Q * 0.2, the negative pressure gradient 
is reduced to a very small value. For m o = 0.4 and higher, the pressure 
gradient becomes positive. Thus, the effect of the surface blowing on a 
sphere-cone is to smooth out the curvature discontinuity at the juncture 
point which results in a positive pressure gradient over the body. 

Shock standoff distances are shown in figure 9. This figure also shows 
the smoothing out of the curvature discontinuity with increasing blowing. 

It is also seen that the shock gets steeper and steeper as the blowing rate 
is increased. It did not create any problem for the present conditions and 
body shape, but the flow can become subsonic behind the shock for large angle 
bodies at downstream points. 

Figures 10 through 12 show the profiles of tangential velocity enthalpy 
and normal velocity at s = 0.19635 for the 45-degree half-angle sphere-cone 
with m Q increasing from 0 to 0.6. Profiles at s - 0.19635 are also plotted 
in figures 13 through 15 for a 45-degree half-angle hyperboloid using the 
second set of conditions with i!: = 0 and 0.4. It is seen that, in the nose 
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region, the profiles for the sphere-cone are very similar to that for the 
hyperboloid. The normal velocity profile for the sphere-cone also shows the 
same behavior as was seen for the hyperboloid, l.e., it decreases slightly In 
the inner Invlscid layer, increases in the thin shear layer, and then decreases 
monotonically to the shock value in the outer Inviscid layer. 

Tangential velocity, enthalpy and normal velocity profiles are plotted 
at s * 1.76715 in figures 16 through 21 for both the sphere-cone arid hyper- 
boloids with increasing m Q . It is seen from figure 17 that, for the sphere- 
cone, the enthalpy changes rapidly in the inner portion of the outer invlscid 
layer forming an entropy layer. For the hyperboloid, the flow quantities 
change gradually to the shock value in the outer inviscid layer. This 
difference in the profiles for the two bodies will have an effect on the 
radiative heating. 

The effect of turbulence on the flow field over the sphere-cone is 
discussed In the presence of surface injection. For the turbulent solutions, 
the flow is assumed to undergo instantaneous transition at s = 0.19635. 

Figure 22 shows the surface pressure distribution for the laminar and turbulent 
flows. For zero blowing, the surface pressures are almost identical in the 
laminar and turbulent flows. For m Q a 0.4, the laminar flow shows a positive 
pressure gradient, whereas the turbulent flow develops again a negative pressure 
gradient similar to that for zero blowing. 

Figure 23 shows the shock standoff distances for laminar and turbulent 
flows. Here again, the shock standoff distances for laminar and turbulent 
flows are almost identical for zero blowing. For nonzero blowing, there is a 
significant decrease in the shock standoff distances for the turbulent flow 
at downstream points. Moreover, the shock slope *or the turbulent flow is 
not seen to increase with increasing blowing rates. The figure shows that, 
for the turbulent flow, the shock is just pushed away from the body in the 
presence of blowing but the shock shape remains similar to that for zero 
blowing. The reduction in the shock standoff distance can be explained 
from figure 16 in which the turbulent tangential velocity profile for m Q = 

0.4 is plotted along with the laminar profiles. The magnitude of the 
tangential velocity is higher all across the shock layer for the turbulent 
flow, which results in the reduction of the shock standoff distances. The 
turbulent tangential velocity profile also shows the formation of a boundary- 
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layer-11 ke region near the surface. Corresponding changes In the turbulent 
enthalpy profile are seen In figure 17 where the enthalpy near the surface 
goes up for the turbulent flow. The turbulent enthalpy profile also shows a 
reduction In the gradients in the entropy layer due to the fact that the 
shock Is relatively less strong for the turbulent flow than for the laminar 
flow. The increased gradients near the surface in the turbulent flow make 
the convective heating and skin friction significant even with massive 
blowing whereas these are negligibly small for the corresponding laminar 
flow. 

Another very important effect of the formation of the boundary-1 ayer- 
1 1 ke region is on the radiative heat flux toward the body. It Is discussed 
in reference 2 for the coupled flow-field solutions over a hyperboloid for 
Jovian entry. Due to the formation of a high-temperature ablation layer 
near the surface, the radiation blocked in the ablation layer is greatly 
reduced for the turbulent flow, which significantly increases the radiative 
heat transfer to the body surface. 

Figure 24 shows the eddy viscosity distribution across the shock layer 
at s = 1.76715 for the 45-degree half-angle sphere-cone with m Q = 0, 0.2, 
and 0.4. It is seen from this figure that the eddy viscosity increases 
sharply with the increasing surface blowing. The two-layer model provides 
good approximation of the eddy viscosity for zero or small blowing rates, but 
the validity of this model with massive surface blowing is not known. 

However, it is shown in reference 2 that even if the two-layer model is 
assumed to be in error by a factor of 10, the radiation heat flux to the 
Jovian entry probe is substantially increased. 

6. CONCLUSIONS 

Numerical solutions for the laminar and turbulent flow over a spherically 
blunted cone and hyperboloid are presented with massive surface blowing. It 
is found that the effect of the surface blowing on the laminar flow field over 
a sphere-cone is to smooth out the curvature discontinuity at the juncture 
point, which results in a positive pressure gradient over the body. The shock 
slope increases at the downstream points with increasing surface blowing. 

The turbulent flow with surface blowing is found to redevelop a boundary- 
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layer-like region near the surface. The Increased gradients near the surface 
In this boundary-layer-1 Ike region make the convective heating and skin 
friction significant even with massive blowing. The shock standoff distances 
are significantly reduced over the downstream part of the body which make the 
shock relatively less strong for the turbulent flow than for the laminar 
flow. 

Results presented here are for the perfect gas-type flow and do not 
Include radiative heating or chemical reactions In the flow field. The 
main purpose of this report Is to demonstrate the capability of the present 
computer code In providing the solutions for the flow over a realistic probe 
shape with massive surface blowing. The code will be modified to provide 
the coupled solutions for the chemically reacting and radiating flow with 
surface ablation. 
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Figure 1. Coordinate system. 
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Figure 13. Tangential velocity profile variation with surface 
blowina at s = 0. 19635. 
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Figure 14. 
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Figure 15. Normal velocity profile variation with surface blowing at 
s = 0.19635. 
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Figure 21. Normal velocity profile variation with surface 
blowing at s = 1.76715. 
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Figure 23. Shock standoff distances for laminar and turbulent 
flows with surface blowing. 
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Figure 24. Eddy viscosity profile variation with surface blowing at 
s = 1.76715. 
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